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ABSTRACT 

Several  anti-aliasing  strategies  are  proposed, 
which  generate  Monte  Carlo  discretized  estimates 
of  color  and  intensity  at  each  pixel  of  a  raster 
display. 
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1.  We  are  given  a  function  f:  R  ■»  R  ,  specifying  color  and  intensity  at 

2 

any  point  of  a  screen  area  S  £  R  .  The  screen  S  is  subdivided  into  .7  r 
?•  (h  =  1,  2,  ....  .7),  all  disjoint  and  of  equal  area  and  shape. 

2.  It  is  intended  to  approximate  the  function  f  on  5  by  a  function  i :  R“  - 

which  takes  the  value  a-  on  the  pixel  ?,  ,  for  h  =  1.  2,  ....  .7. 

n  r  h 

5.  One  approach  is  to  define,  for  the  pixel  centered  at  c :  a  u'eig'n r 

-unction  \'{r  -  c.  )  =  -j,  ( r)  and  let 
k  n  h 


r 

^  3  |  dr  f(r)uh(r), 
'  2 


a) 


where  £  denotes  R“  and  f„dr  denotes  f°°  d xf°°  du ,  with  r  =  (x,  >). 

^  ^  '  -  CO  '  —00  w'  *■ 

4.  A  very  general  Monte  Carlo  scheme  for  estimating  ?.  would  select 

n 

integer  n.  and  a  set  of  estimator-probability  pairs  (0*  o-  .  (  O  )  » 

f 1  nZ  ^Z 

:  =  1,  2,  ....  n,;  so  that  one  samples  points  f  .  6  £  with  probability 
independently  of  each-other,  and  uses  the  csCimal^r 
nh 


l  9hA$;) 

•  |  '  L  is  is 


(2) 


an 

for 

dens  i  t 


for  For  example,  "crude  Monte  Carlo"  could  define  o-  • ( f)  =  .7/4 ,  where 

A  is  the  area  of  2  (so  that  A/N  is  the  area  of  the  pixel  P.) ,  and  use  the 
estimator  g^(r)  =  af{r )  in  P.^\  but  this  would  not  work,  since  we  would 
want  that  the  estimator  be  unbiased,  i.°.,  that 


t  =  l  ^  m  n 

and  this  reduces,  by  f  1  j  , '  to  c  =  At^/Nb^n^,  w'.,cre 


(3) 


d'  f{r)  , 


and  we  would  need  to  know  both  t-,  and  to  get  j.  !  Anotner  approach  is 

to  use  a  -■{.*)  -  uv(r)  =  a'  ( r  -  cj  in  the  whole  of  i  (.though,  of  course, 
most  of  the  probability  will  be  in  or  near  t^)  ,  and  use  the  estimator  ) 

-  jf( r) ;  whence  the  condition  (3)  reduces  to  2  =  1/K>t  provided  that  the 
weight  function  satisfies  (.as  is  usual)  the  normalicing  condition 


d r  a.  (  r)  =  drj(r-  c,  )  =  d r  u (  r)  =  1 , 


Of  course,  this  condition  is  not  at  all  unreasonable.  Note  that  we  may, 

yet  again,  choose,  over  the  whole  of  Q,  o._.(r)  =  ( r) »  a  different  normalized 

fl'l, 

weight  function  from  (for  instance,  the  normal  distribution  centered  ate. 
and  with  standard  deviation  of  the  order  of  the  diameter  of  a  pixel),  and 
then  the  estimator  would  be  g^.{r)  =  r)f{  r)/u^(  r)n^,  as  is  readily 
verified,  and  this  is  again  feasible;  so  we  note  the  pair: 

U(0/(r)  ) 

^9hi ”  =  1  u/C r)n,  *  w  '6-) 

n  n 

3.  An  alternative  approach  would  be  to  use  a  form  of  szvzzif-sd  our: r’ivp. 
Note  that, in  the  technique  developed  above,  all  n ,  estimators  are  identical 
and  identically  distributed.  Suppose,  instead,  that  the  pixel  P,  is  dis¬ 
sected  into  m  identical  sub-pixels  /?,  .,  and  that  s  .  identical  estimators 

K  hj  J 

y,  .ft)  nr°  sampled  with  density  p,.(r)  in  Q,  where  .  ( r)  -  o(r  -  h  )  and 

ri,j  rlQ  KlQ 

b.  .  is  the  center  of  R.  ..  We  then  require,  by  (3),  that 


-3- 


l  =  .  d r  f(r)u  (r) . 


,-  =  l  v  JQ 


As  an  example,  we  could  choose  the  function  o,  and  then  put 


f(r)  v(r  -  C) 


^  r ^  ms  .  c  (  r  -  6,  .)  ' 


where  we  also  must  have  that 


s  .  =  n. 


6.  IvTiat  we  must  do  to  make  the  method  efficient  is  to  minimize  (or  at  least 
diminish)  the  variance  of  our  estimate.  Thus,  we  note  that,  for  the  first 
technique,  given  by  (6),  we  have 


i  \~h v  v 

££i  vari3«]  ■  nnQir[-^h 


—  rrr; -  U,'  (  r 


w^nnh  *h 


Wr)  }mh\‘*h2)’  (10) 

1  1  n 


where 


h^inrifinv 


For  the  second  technique,  given  by  (8),  we  similarly  get  that 

.  ~  i  ff  (/('M'  -  c  ),2 

var>  >  s  .g,  .}  =  >  e.  var [g.  .]  =  V  s  .<  dr - - - r— —  o{r  -  b  ) 

Ai  J  ^  A  •'ho  K-» ( '  -  A ,o I  'l 


f(  ')ul  r  -  c, ) 

9dr=?‘'-  V"('*  V 


X  “T- f“itj  '  ',(i2)  ■ 

j=\  m  8  .  ° 

J 
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where 


..*(/■)  ]“>(  r  -  c.  ) 

_dr  p(r  -  b..) 


(13) 


If  we  consider  the  case  of  (6),  (10),  and  (11),  and  first  assume  that 
.* »  '  >  and  so  2>;.  and  are  all  given  a  priori',  then  we  may  ask  how  to 

choose  the  numbers  of  function-evaluations  n,  by  pixels,  so  as  to  make  all 
variances  the  same,  given  the  sum  n  =  lN.  n.  .  The  answer  is  evidently 

nh  -  -  **2>.  (i-D 

and  the  common  value  of  the  variance  at  every  pixel  is  then 

vari-'d  =  l  -  =  ^  )/n-  (15) 

In  the  case  of  (8),  (12),  and  (13),  with  f ,  u ,  p,  and  so  and  given,  we 
similarly  see  that  we  can  first  optimize  over  the  strata  in  a  single  pixel; 
Lagrangian  theory  shows  that 


s* =  Mv*  -  ♦.  -ifru,  -  h2)u 


minimizes  the  variance  at  P.  to  the  value 

n 


min  var'L  l  .  ,  a,  .1  = 

-  =  \  ~  r*  *  *  2  I 

c  n  n*  ^  ° 

h. 


1  [jn  ,  .  2^1 2 

2.  rj*iCuW  '  *k  ]  * 


(16) 


cm 


Note  that  the  Cauchy-Schwart z-Bunyakovsky  inequality  shows  that  indeed 


1  v’77 


1  Lm  ,  2  2  _ 

2  rj  =  iuhj  J  j  '  2  V=i 

n  h 


tu  -y.1, .  n 
1  • 


s  . 


J _ (zw 


h  )  E? 


2„  i  J  =  1  8.  j  1=1  1 


(18) 


m  n 
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and  the  right-hand  side  of  the  inequality  is  the  general  variance  (12),  by 
(9);  so  that  (16)  does  indeed  minimize  (not  maximize  or  point-of-inflexion) 
the  variance.  Now  we  proceed,  as  before,  to  make  all  the  variances  (17)  the 


samQ;  yielding  that 


-  !  '  2  r 


- >  i , :  “ 


-  ,  IP'  •  * 

c  !  ^=1  nj  n 


KUh.i  ‘  J  /Lk=l\L-'=lV'ik,;  " 


This  makes  the  common  value  of  the  variance 


min  var[Z 


1  JV  L  m 


J'=l  ~  m2n 


2^K 


^il^-i^  -  vn  • 


8.  As  a  specific  example,  we  may  suppose  that  5  is  a  rectangle 

5  =  (0  <  x  <  Lv  0  <  y  ^  L2);  (21) 

and  that  the  index  h  is  (^  ,  with  N  =  .V^  and  0  *£  h.^  <  ,V+  (t  =  1,  2), 

so  that  ?,  is  the  rectangle 


P.  =  ?, 
n  h 


^  --  +  +  f22> 


centered  at  c  =  (okr  3jj,)  with  =  rf-  (&*  ♦  h)  (t  =  1,  2). 


Similarly,  we  take  j  =  (j^ ,  j*2) ,  m  =  ny«2,  and  0  (t  =  1,  2),  so 

that  ..  is  the  )  rectangle  centered  at 

L 

=  (b,  ,  o,  ._)  with  b,  . ,  =  — — - — (m  b  +7  «  Ui  f  +■  =  i  T\  (7*i 
'V  -vi’  bj2  bji  ;;  w  1  t  6  Jt  }  1 

We  may  further  postulate  that  both  and  p  take  the  form  of  the  normal 
distribution,  with 


-6- 


I  r‘  =  3^7  exp(-{(x  *  ^;;1)2  *  (y  -  ^2)2}/2y),  (25- 

where  y  =  (^I,/.'.'^)  a  =  («/.V)c,  (26) 

and  z;,..(r)  =  exp(-'(x  -  ru.,.,)2  +  (y  -  b.  ..,)  ^;/2£)  .  (2‘) 

't/  *• -l  "v  ■*  ,4,t/  ^ 

where  E  =  (A/’ln  r.^)a  =  (4/.7m)o.  (28) 

Here,  :  is  a  constant  for  the  system,  related  to  the  weight  function  u  but 
not  to  f  or  to  5  and  its  subdivisions. 

Then  we  have  that 


I  i  i~">  •)  , 

dx ,  dy  *’(x,  u))  <-j(x  -  a,  ,  ,  v  -  c?,  ,) 

0  Jo  “  *  rtl  “  ' 

x  exp ( j{  (x  -  ahl)2  *  (y  -  o^ri/Ea) 


(29) 


anci  ^  2-a  '  1  dxj  2  d y  [f(x,  y))2[w{x  -  j  -  e.,.,)]2 

17  •  0  J  0 


>,1‘ 


>.2 


x  exp (^{  (x  -  b-  ..)2  +  (y  -  £.  „.,) 2 > / 2 a )  .  (30) 

9.  The  strategies  investigated  here  so  far  are  adaptive  only  insofar  as 

the  optimizing  numbers  of  samples  (14)  and  (16)  are  to  be  estimated  from 

Monte  Carlo  estimates  of  the  and  u,  .which  can  be  obtained  s imul taneous lv 

n  nj 

with  the  estimates  of  ^  generated  by  the  estimators  (6)  and  (8),  respectively. 
Since  only  small  samples  are  to  be  taken,  because  f  is  so  laborious  to  get, 
the  relative  sample-sizes  (14)  and  (16)  will  not  be  very  accurately  optimal. 


Another  approach  would  attempt  to  perform  importance  sampling  by 
sequentially  approximating  f(x,  y)w^(x,  a)  with  w'.  Sincec\is  given  and 
f  is  experimentally  determined  (so,  also  given),  we  may  write  T(x,  y)  for 
the  product.  As  we  accumulate  values  of  C  by  sampling  (initially  with  an 
arbitrary  distribution),  we  can  form  an  increasingly  accurate  picture  of 
the  functional  dependence  of  C  on  (x,  y)  and  model  y/  on  this. 

Alternatively,  we  may  do  a  sequential  correlated  sampling  calculation, 
in  which  we  fix  the  sampling  density  arbitrarily,  and  then  use  an  estimator 
of  the  form  (C(x,  u )  -  i/(x,  y)}/w.'(x,  >)  -  fndr  y(x,  u) ,  where  is  the 
best  approximation  to  f  for  which  the  integral  on  the  right  is  easily 
computable, 

10.  Yet  another  approach  which  should  be  empirically  investigated  is  to 
use  an  ordering  of  the  sampled  values  of  C  to  indicate  where  stratification 
should  occur.  First,  we  sample  C  at  a  small  number  of  points  in  each  pixel 
and  tabulate  C,  x,  and  y,  in  order  of  increasing  C.  If  there  is  a  strong 
correlation  of  C  with  x  or  with  y,  split  the  pixel  accordingly  and  sample 
a  few  more  points.  Repeat,  if  necessary, 

dote  that  the  strati fication  and  sampling  are  dene  in  the  wo tie 
not  wi  thin  the  pixel  or  sub-pixel  only.  This  is  to  conform  with  the  global 
for m  of  w.  dote  also  that  w  may  be  given  the  full  theoretical  form,  and 
need  not  be  approximated  by  a  normal  distribution  itself. 


